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Recently, important updates were made for the hadronic contribution to the theoretical pre- 
diction of g — 2. The isospin-breaking-corrections, needed in the comparison of the two pion 
spectral functions from r decays and e'^e~ annihilations, were improved using new exper- 
imental and theoretical input. The recently published BABAR data were included in the 
global average of e'^e~ spectral functions. These data, as well as the ones from r decays, were 
combined using newly developed software, featuring improved data interpolation and aver- 
aging, more accurate error propagation and systematic validation. The discrepancy between 
the e'^e~ and the r-based result is reduced from previously 2.4 to 1.5 a. The full Standard 
Model prediction of g — 2, obtained using e'^e~ data, differs from the experimental value by 
3.2 standard deviations. 



1 Introduction 

The Standard Model (SM) prediction for the anomalous moment of the muon can be conviniently 
separated into a sum, 

= "^''^ + + «r'' (1) 

with the QED, hadronic and weak contributions respectively. This prediction is limited in 
precision by the lowest-order hadronic vacuum polarization contribution, which together with 
the hadronic higher order and light-by-light (LBL) contributions provide the total hadronic 
contribution. Owing to unitarity and to the analyticity of the two point correlator, using the 
optical theorem, the lowest order hadronic vacuum polarisation contribution to the anomalous 
magnetic moment of the muon {a^^'^'^'~') can be computed through an energy-squared dispersion 
integral (ranging from the n^j threshold to infinity): 
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where K{s) is a QED kernel function [I]. Actually, the integration kernel strongly emphasises 
the low-energy part of the spectrum, about 73% of the lowest order hadronic contribution being 
provided by the 7r7r(7) final state §. More importantly, 62% of its total error stems from the 
7r7r(7) mode, stressing the need for ever more precise experimental data in this channel to confirm 
or not the observed deviation between SM prediction and experiment. 

In this note we concentrate on the computation of a^^'^'^'~' . We add the other contributions 
and compare the total SM prediction with the experimental measurement. The results presented 
here were first published in Ref. [2J. We also recall some previous results from Ref. [3]. 



" Throughout this note, final state photon radiation is implied for hadronic final states. 



2 Hadronic data and contributions to ' 

A former lack of precision e^e~ -annihilation data inspired the search for an alternative. It 
was found [4J in form of r — )• z^,- + 'k~tt^ , 27r^7r+7r'^, tt^Stt'^ spectral functions [SHS], transferred 
from the charged to the neutral state using isospin symmetry. During the last decade, new 
measurements of the tt~^tt~ spectral function in e^e~ annihilation with percent accuracy became 
available [9Hl2] , superseding or complementing older and less precise data. 

Indeed, fix-energy measurements from the CMD2 pTOJ and SND [llj experiments at the 
VEPP-2M collider (Novosibirsk, Russia), achieved comparable statistical errors, and energy- 
dependent systematic uncertainties down to 0.8% and 1.3%, respectively. These measurements 
have been complemented by results from KLOE [12j, where a hard-photon ISR technique was 
applied for the first time to precisely determine the tt+tt^ cross section between 0.592 and 
0.975 GeV. The analysed data sample provides a 0.2% relative statistical error on the vr+vr" 
contribution to a^'^^ . KLOE does not normalise the 7r"'"7r~7 cross section to e^e'~' — t- /i^^~7 
so that the ISR radiator function must be taken from Monte Carlo simulation (c/. 1.13j and 
references therein). The systematic error assigned to this correction varies between 0.5% and 
0.9% (closer to the (/> peak). The total assigned systematic error lies between 0.8% and 1.2%. 

With the increasing precision, which today is on a level with the r data in that channel, 
systematic discrepancies in shape and normalisation of the spectral functions were observed 
between the two systems |141I15|. It was found that, when computing the hadronic VP con- 
tribution to the muon magnetic anomaly using the r instead of the e^e~ data for the 27r and 
47r channels, the observed deviation with the experimental value [T6] would reduce to less than 
la ^17j. The discrepancy between the r and e"'"e~-based predictions decreased after the inclu- 
sion of new T data from the Belle experiment |18] . published e^e~ data from CMD2 |10j and 
KLOE [12] (superseding earlier data Jl9] ) . and a reevaluation of isospin-breaking corrections 
affecting the r-based evaluation [3]. Q In terms of a^'^'^'^ , the difference between the r and 
e^e~-based evaluations in the dominant tt+tt" channel was found to be 11.7ib3.5ee ±3.5t-+ib [3] 
(if not otherwise stated, the values are given in units of lO"^*'), where KLOE exhibits the 
strongest discrepancy with the r data (without the KLOE data the discrepancy reduces from 
2.4cr to 1.9<t). 

Not long time ago the BABAR Collaboration reported [20j measurements of the processes 
e+e~ — )• 7r+7r~(7), /i"'";^~(7) using the ISR method at 10.6 GeV centre-of-mass energy. The 
detection of the hard ISR photon allows BABAR to cover a large energy range from threshold 
up to 3 GeV for the two processes. The 7r+7r~(7) cross section is obtained from the -7r"'"7r~7(7) to 
/x"'"/i~7(7) ratio, so that the ISR radiation function cancels, as well as additional ISR radiative 
effects. Since FSR photons are also detected, there is no additional uncertainty from radiative 
corrections at NLO level. Experimental systematic uncertainties are kept to 0.5% in the p peak 
region (0.6-0.9 GeV), increasing to 1% outside. 

3 Combining cross section data 

The requirements for averaging and integrating cross section data are: {i) properly propagate all 
the uncertainties in the data to the final integral error, [ii) minimise biases, i.e., reproduce the 
true integral as closely as possible in average and measure the remaining systematic error, and 
{in) optimise the integral error after averaging while respecting the two previous requirements. 
The first item practically requires the use of pseudo-Monte Carlo (MC) simulation, which needs 
to be a faithful representation of the measurement ensemble (with all known correlations) and to 
contain the full data treatment chain (interpolation, averaging, integration). The second item 

^ The total size of the isospin-breaking correction to aj^^'*''"'"' has been estimated to (— 16.1 ± 1.9) • 10"^" , which 
is dominated by the short-distance contribution of (—12.2 ± 0.2) • 10~^° [3]- 



requires a flexible data interpolation method (the trapezoidal rule is not sufficient as shown 
below) and a realistic truth model used to test the accuracy of the integral computation with 
pseudo-MC experiments. Finally, the third item requires optimal data averaging taking into 
account all known correlations to minimise the spread in the integral measured from the pseudo- 
MC sample. The details of our newly developed combination procedure (HVPTools) are given 
in Ref. [2]. 

HVPTools transforms the bare cross section data and associated statistical and systematic 
covariance matrices into fine-grained energy bins, taking into account to our best knowledge the 
correlations within each experiment as well as between the experiments (such as uncertainties in 
radiative corrections). The covariance matrices are obtained by assuming common systematic 
error sources to be fully correlated. To these matrices are added statistical covariances, present 
for example in binned measurements as provided by KLOE, BABAR or the r data, which are 
subject to bin-to-bin migration that has been unfolded by the experiments, thus introducing 
correlations. 

The interpolation between adjacent measurements of a given experiment uses second order 
polynomials. This is an improvement with respect to the previously applied trapezoidal rule, 
corresponding to a linear interpolation, which leads to systematic biases in the integral. In the 
case of binned data, the interpolation function within a bin is renormalised to keep the integral 
in that bin invariant after the interpolation. The final interpolation function per experiment 
within its applicable energy domain is discretised into small (1 MeV) bins for the purpose of 
averaging and numerical integration. 

The average weights of the interpolated measurements from different experiments contribut- 
ing to a given energy bin takes into account the precision and the correlations of different 
measurements and experiments, as well as different measurement densities or bin widths within 
a given energy interval. This provides an optimal final uncertainty on a^^'^'^'^ . 

The left hand plot of Fig. [1] shows the weights versus ^/s the different experiments carry in 
the average. BABAR and KLOE dominate over the entire energy range, completely covered by 
BABAR. 

If the value of a bin- wise average exceeds the number of degrees of freedom (n^of ) , the error 
in this averaged bin is rescaled by \/ X^/^-dof to account for inconsistencies (c/. Fig. [H Right). 
Such inconsistencies frequently occur because most experiments are dominated by systematic 
uncertainties, which are difficult to estimate. 

We have tested the fidelity of the full analysis chain (polynomial interpolation, averaging, 
integration) by using as truth representation a Gounaris-Sakurai vector-meson resonance model 
faithfully describing the vr+vr" data. The difference between true and estimated a^^'^'^'^ values 
is a measure for the systematic uncertainty due to the data treatment. We find negligible bias 
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Figure 1: 
the figure 



Left: relative averaging weights per experiment versus ^/s. The experiments labelled "other exp" in 
correspond to older data with incomplete radiative corrections. Right: rescaling factor accounting for 
inconsistencies among experiments versus ^/s (see text). 
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Figure 2: Relative cross section comparison between individual experiments (symbols) and the HVPTools average 
(shaded band) computed from all measurements considered. Shown are BABAR (left) and KLOE (right). 



Table 1: Evaluated aj^'^'*''^ [tttt] contributions from the e^e data for different energy intervals and experiments. 
Where two errors are given, the first is statistical and the second systematic. We also recall the r-based result. 



Energy range (GeV) 


Experiment 


aJl^'l'LOM (10-10) 


2m^± - 0.3 


Combined e+e" (fit) 


0.55 ±0.01 


0.30 - 0.63 


Combined e+e" 


132.6 ± 0.8 ± 1.0 (l.Stot) 


0.63 - 0.958 


CMD2 03 


361.8 ±2.4 ±2.1 (3.2tot) 




CMD2 06 


360.2 ±1.8 ±2.8 (3.3tot) 




SND 06 


360.7 ±1.4 ±4.7 (4.9tot) 




KLOE 08 


356.8 ±0.4 ±3.1 (3.1tot) 




BABAR 09 


365.2 ±1.9 ±1.9 (2.7tot) 




Combined e~^e~ 


360.8 ±0.9 ±1.8 (2.0tot) 


0.958 - 1.8 


Combined e~^e~ 


14.4 ±0.1 ±0.1 (0.2tot) 


Total 


Combined e^e~ 


508.4 ±1.3 ±2.6 (2.9tot) 


Total 


Combined t 


515.2 ± 2.0cxp ± 2.28 ± l-9iB (3.5tot) 



below 0.1 (remember the 10"^'' unit), increasing to 0.5 (1.2 without the high-density BABAR 
data) when using the trapezoidal rule for interpolation instead of second order polynomials. 

Fig. [2] shows the relative differences between BABAR, KLOE, and the average. Fair agree- 
ment is observed, though with a tendency to larger (smaller) cross sections above ~0.8 GeV for 
BABAR (KLOE). It is interesting to point out that the same type of tendency shows up when 
comparing the KLOE measurement with the IB-corrected r average [3J. These inconsistencies 
(among others) lead to the error rescaling shown versus ^/s in Fig. [H Right. 

A compilation of results for a|l^'^'^^ [vrvr] for the various sets of experiments and energy regions 
is given in Table [TJ The comparison with our previous result [3], aji^'^'^*^ [tttt] = 503.5 ± 3.5tot5 
shows that the inclusion of the new BABAR data significantly increases the central value of the 
integral, without however providing a large error reduction. This is due to the incompatibility 
between mainly BABAR and KLOE, causing an increase of the combined error. In the energy 
interval between 0.63 and 0.958 GeV, the discrepancy between the aji*^*^'^^ [tttt] evaluations from 
KLOE and BABAR amounts to 2.0fT. Using only the BABAR data to evaluate aji'^'^''^*^[7r7r] one 
finds [20J 514.1 ± 2. 

2stat i 3.1syst, which is in very good agreement with the result from the t 
average. Including BABAR in the global average reduced the difference between the r and e^e~- 
based predictions to l.ba. A difference (slope) in the shape of the two spectral functions subsists 
however [2], mainly due to the difference between the r average and KLOE measurement. 

We also reevaluate the e^e~ — )• ■n~^i:~2TT^ contribution to a^'^^ . The CMD2 data used 
previously [21j have been superseded by modified or more recent, but yet unpublished data [22], 
recovering agreement with the published SND cross sections [23]. Since the new data are unavail- 



able, we discard the obsolete CMD2 data from the vr+vr' average, finding ' [7r7r27r''] = 
17.6 ±0.4 

Stat i l-7syst (compared to 17.0 i 0.4stat i l-6syst when including the obsolete CMD2 

data). 

Adding to the e"'"e~-based Uff ' [vrvr] and ' [7r7r27r'^] results the remaining exclusive 
multi-hadron channels as well as perturbative QCD |17j . we find for the complete lowest order 
hadronic term a^^'^'^ [e^e~] = 695.5 it 4.0cxp ± 0.7qcd (4.1tot) • It is noticeable that the error 
from the tt'^tt" channel now equals the one from all other contributions to a^^'^'^'^ . Adding 
further the contributions from higher order hadronic loops, hadronic light-by- light scattering, 
as well as QED and electroweak effects, we obtain the SM prediction (still in 10"^*^ units) 

af^[e+e-] = 11 659 183.4 ± 4.1 ± 2.6 ± 0.2 (4.9tot) , 

where the errors have been split into lowest and higher order hadronic, and other contributions, 
respectively. The a^'^[e'''e~] value deviates from the experimental average by 25.5 it 8.0 {3.2a). 
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Figure 3: Compilation of recent results for a^, subtracted by the central value of the experimental average |25] . 
The shaded vertical band indicates the experimental error. The SM predictions are taken from: HMNT 07 |24) . 
JN 09 [26], Davier et al. 09 3 (r-based and e^e~ including KLOE), and the e^e~-based value from this work. 



A compilation of recent SM predictions for compared with the experimental result is 
given in Fig. [3l The BABAR results are not yet contained in evaluations preceding the present 
one. The result by HMNT ^24j contains older KLOE data |19J, which have been superseded by 
more recent results [T2], leading to a slightly larger value for aji^*^'^^. 



4 Conclusions 



We have presented a reevaluation of the lowest order hadronic contribution to the muon magnetic 
anomaly in the dominant vr+vr^ channel, using new precision data published by the BABAR 
Collaboration. After combination with the other e~^e~ data a 1.5cr difference with the r data re- 
mains for the dominant vr^vr" contribution. For the full e"'"e~-based Standard Model prediction, 
including also a reevaluated 7r^TT^2'K^ contribution, we find a deviation of 3.2a from experiment 
(reduced from 3.7a without BABAR). The deviation reduces to 2.9a when excluding KLOE 
data, and further decreases to 2.Aa when using only the BABAR data in the vr+vr" channel. As 
a reminder, the r-based result deviates by 1.9a from the Standard Model. 

The present situation for the evaluation of aji'^'^''"'^[7r7r] is improved compared to that of recent 
years, as more input data from quite different experimental facilities and conditions have become 
available: e^e~ energy scan, e"'"e~ ISR from low and high energies, r decays. Our attitude has 
been to combine all the data and include in the uncertainty the effects from differences in the 
spectra. At the moment the ideal accuracy cannot be reached as a consequence of the existing 
discrepancies due to uncorrected or unaccounted systematic effects in the data. A critical look 



must be given to the different analyses in order to identify their weak points and to improve on 
them or to assign larger systematic errors. 

Problems also persist in the 7r"^7r~27r*^ mode, where the r and e^e~-based evaluations differ 
by (3.8±2.2)-10-^°, but also the e+e" data among themselves exhibit discrepancies. Fortunately, 
new precision data from BABAR should soon help to clarify the situation in that channel. 
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